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ABSTRACT 

We discuss how the space of possible cosmological parameters is constrained by the angular diameter 
distance function, Da{z), as measured using the SZ/X-ray method which combines Sunyaev-Zel'dovich 
(SZ) effect and X-ray brightness data for clusters of galaxies. New X-ray satellites, and ground-based 
interferometers dedicated to SZ observations, should soon lead to Da(z) measurements limited by sys- 
tematic rather than random error. We analyze the systematic and random error budgets to make a 
realistic estimate of the accuracy achievable in the determination of (il m , ^A, h), the density parameters 
of matter and cosmological constant, and the dimensionless Hubble constant, using Da(z) derived from 
the SZ /X-ray method, and the position of the first "Doppler" peak in the cosmic microwave background 
fluctuations. We briefly study the effect of systematic errors. We find that f2 m , f^A, and w are affected, 
but h is not by systematic errors which grow with redshift. With as few as 70 clusters, each providing a 
measurement of Da{z) with a 7% random and 5% systematic error, Q m can be constrained to ±0.2, Q\ 
to ±0.2, and h to ±0.11 (all at 3a). We also estimate constraints for the alternative three-parameter set 
(Cl m , w, h), where w is the equation of state parameter. The measurement of Da{z) provides constraints 
complementary to those from the number density of clusters in redshift space. A sample of 70 clusters 
{Da measured with the same accuracy as before) combined with cluster evolution results (or a known 
matter density), can constrain w within ±0.45 (at 3tr). Studies of X-ray and SZ properties of clusters of 
galaxies promise an independent and powerful test for cosmological parameters. 



1. INTRODUCTION 

What set of cosmological parameters characterizes our 
Universe? 

According to the most popular cold dark matter (CDM) 
scenario, the Universe consists of baryonic matter and a 
substantial amount of "dark" matter. A variety of recent 
measurements have led to the conclusion that the matter 
density parameter fl m ss 0.3 (Turner 2000), while CMB 
measurements strongly favor a flat space time (Bernardis 
et al. 2001; Lee et al. 2001), and SNe la measure- 
ments indicate that the Universe is accelerating, suggest- 
ing a negative pressure (Riess et al. 2000; Perlmutter 
et al. 1999). Taken together, these pieces of evidence 
suggest that the baryonic and dark matter content of the 
Universe is supplemented by an additional smooth compo- 
nent with negative pressure, P w , modeled by the equation 
of state p w c 2 = —wP w , where p w is the density of this 
component, w is a dimensionless state parameter of order 
unity (cf. Huterer & Turner 2000). 

Each existing dataset constrains, with limited accuracy, 
some subset of the cosmological parameters. Different 
measurements and combinations of measurements, such 
as SNe la, Cosmic Microwave Background (CMB) fluctua- 
tions, IRAS infrared galaxy surveys, classical double radio 
galaxy properties, 1.2-Jy galaxy redshift surveys, gravi- 
tational lensing, cluster X-ray temperature function and 
cluster number counts, baryon and gas mass fraction, and 
the SZ effect have been used to constrain cosmological pa- 



rameters (Jaffe et al. 2000; Balbi et al. 2000; Tegmark & 
Zaldarriaga 2000; Guerra, Daly & Wan 2000; Efstathiou 
et al. 1999; Lasenby et al. 1999; Perlmutter et al. 1999; 
Gawiser & Silk 1998; Lineweaver 1998; White 1998; Pen 
1997; Sasaki 1996; Huterer & Turner 2000; Majumdar 
& Subrahmanyan 2000; Bridle et al. 1999; Diego et al. 
2001). 

In the future, the SNAP project (http://snap.lbl.gov) 
plans to use the SNe la method to determine the matter 
density and the cosmological constant at the few percent 
level. Even with the next-generation CMB satellites, MAP 
and Planck, degeneracies will remain among the cosmo- 
logical parameters that can be estimated from the results 
(Efstathiou & Bond 1999; Zaldarriaga, Spergel & Seljak 
1997). The importance of using a wide range of methods, 
therefore, is twofold. First, a simultaneous consideration 
of all data sets should allow the best joint estimation of the 
cosmological parameters. Second, the agreement of dif- 
ferent techniques for measuring the cosmological param- 
eters should provide a cross-check of our understanding 
of the underlying processes and a control against system- 
atic errors. As we extend our analysis of the CMB to 
more complicated models (tensor fluctuations, finite neu- 
trino masses, etc.) the number of cosmological parameters 
increases, and it becomes even more important that the 
widest possible range of datasets is used, and that strong 
controls against systematic errors are in place. 

Many of the techniques of cosmological parameter esti- 
mation use clusters as tracer particles. As a result there is 
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a large number of planned cluster surveys in the two most 
important non-optical observational indicators of cluster- 
ing: the SZ effect and cluster X-ray emission. Sunyaev- 
Zel'dovich effect surveys with dedicated interferometers 
or receiver arrays will observe hundreds of clusters with 
z > 0.5 per year (Browne et al. 2000; Holder, Carlstrom 
& Mohr 2000; Bartlctt 2000). The new X-ray missions 
(Chandra, XMM) will provide data on hundreds of clus- 
ters with high redshift through their deep and medium- 
deep surveys. 

Cluster evolution, the redshift distribution of clusters 
from SZ and X-ray surveys, Nsz(z) and Nx(z), and clus- 
ter number counts as a function of X-ray flux, Nx(S), are 
important constraints on cosmological parameters. While 
methods based on the CMB power spectrum and SNc 
la are sensitive to the angular diameter distance, cluster 
evolution (and number counts) is sensitive to the growth 
function of matter density fluctuations. Bartlett (2000) 
estimated the performance of ground-based, arcminute- 
resolution, SZ surveys and concluded that more clusters 
will be detected with deep, small-area surveys than shal- 
low, wide-area surveys. Kneissl et al. (2001) studied the 
performance of the Arcminute MicroKclvin Imager exper- 
iment and showed that a set of only about 20 clusters, 
with redshifts in the range z = - 0.8 is needed to measure 
Nsz{z) sufficiently well to distinguish between TO = 1 and 
fl m = 0.3 cosmologies. Carlstrom et al. (2001) discuss 
a deep SZ ground based survey, and quantify constraints 
from Nsz{z) on Q m and Nsz{S) and Nsz(z) were 
estimated from the proposed shallower, but all-sky Planck 
survey by Diego et al. (2001), who concluded that about 
300 clusters (with the necessary optical follow-up to mea- 
sure redshifts) would suffice to distinguish between open 
O m = 0.3 and flat Q m = 1 cosmologies at 3a confidence. 
Holder, Haiman & Mohr (2001) discussed the constraints 
on the parameter space defined by (fl m , ft\, as) (where as 
is the normalization of the matter power spectrum) using 
cluster evolution. Holder et al. showed that constraints 
from cluster evolution and SNe la observations are highly 
complementary to each other. Haiman et al. (2000) dis- 
cussed the constraints on the (fi m , w, h) parameter space, 
assuming a spatially flat geometry (Qa = 1 — fi m ), that 
follow from an SZ effect survey and a large angle deep X- 
ray survey (the Cosmology Explorer; Packer and Lamb). 
They found that N$z(z) and Nx(z), combined with con- 
straints from CMB or SNe la experiments, significantly 
reduce the degeneracies between fi m , w, and h. Huterer 
& Turner (2000) estimated the constraints on Q TO and w 
for flat geometry that can be gained by combining results 
from SNAP, Planck and SZ and X-ray surveys. 

As has been realized, the angular diameter distance- 
rcdshift relation, Da(z), is at the heart of many of these 
techniques, and is sensitive to some important combina- 
tions of cosmological parameters while being degenerate 
under others (Jaffe et al. 2000; Tcgmark & Zaldarriaga 
2000; Efstathiou et al. 1999; Lasenby et al. 1999; Perl- 
mutter et al. 1999; White 1998). Recently White (1998) 
estimated constraints on the pairs of quantities (fi m ,fiA) 
and (fl m , w) (the latter in a flat Universe) from the Da(z) 
function based on current SNe la data combined with 
CMB first peak constraints. The analysis shows that the 
constraint on parameters based on Da(z) is nearly or- 



thogonal to the constraint based on the position of the 
first peak in the CMB fluctuation spectrum. These two 
datasets are thus highly complementary, and form a par- 
ticularly powerful pair of measurements (see also Tcgmark 
et al. 1998). 

The shape and normalization of the observed angular 
diameter distance function constrains several cosmolog- 
ical parameters (the standard formulae for distance in 
Fricdmann-Robertson- Walker Universes are given in, for 
example, Peebles 1993). The distance - redshift func- 
tion, Da(z), in CDM models depends on the matter den- 
sity, cosmological constant and Hubble constant, and any 
other particle density which contributes to the curvature 
of space-time. The slope of the distance-redshift function 
at low redshift is a measure of the Hubble constant, while 
the shape of the function depends on the curvature and 
the different densities. In Figure 1 we show the fractional 
difference in Da(z) with fixed matter density and Hubble 
constant (fi m = 0.3, h = 0.65), but various values cos- 
mological constants (Qa = 0.7, 0.6, 0.3, solid, dashed and 
dash dotted lines) relative to a model with zero cosmolog- 
ical constant (Qa = 0). It can be seen from this figure, 
that Da(z) is most sensitive to the value of the cosmolog- 
ical constant at redshift about unity, and quite insensitive 
to that at small or high redshifts. In the redshift interval 
from z = 0.5 to z — 1.8 the angular diameter distance for 
flat (0 TO = 0.3, Cl\ = 0.7) model is more than 10% differ- 
ent from a model with the same matter density but zero 
cosmological constant (Q m = 0.3, Qa = 0). 

We can expect high precision data from hundreds of 
clusters of galaxies in the near future. With the present 
instrument suite, the statistical errors on individual mea- 
surements will be small, and so the usefulness of the data 
will be limited by their systematic errors. In this paper we 
evaluate the error budget of distance determination based 
on the SZ effect and X-ray measurements (assuming that 
the X-ray output is dominated by thermal bremsstrahlung, 
as is appropriate for the hot clusters in which the SZ ef- 
fect is strong), and provide a realistic estimate of errors 
achievable in the angular diameter distance. We estimate 
how well one will be able to constrain the two parameter 
sets (f2 m , Oa, h), and (O m , w, h) (assuming a spatially flat 
geometry, 17a = 1 — £l m ) if this Da{z) function is com- 
bined with the position of the first Doppler peak in the 
angular power spectrum of CMB fluctuations. 

Our treatment is complementary to previous work of 
Kneissl et al. (2001), Carlstrom et al. (2001), Diego et 
al. (2001), Holder et al. (2001), Haiman et al. (2000), 
and Huterer & Turner (2000), who used N sz (z), N SZ (S), 
and Nx(z) to constrain cosmological parameters, since we 
discuss the importance of Da(z). It also complements the 
work by White (1998), who used Da(z) determined from 
existing SNe la data to constrain (Q, m ,0,\), and (Q m ,w) 
(flat): the errors from the SZ/X-ray technique have sig- 
nificantly different characteristics, and we are concerned 
with the limitations that will be encountered with future 
survey data. 

In the next section we briefly describe the well-known 
method of angular distance determination based on the 
SZ effect and thermal bremsstrahlung with an emphasis 
on how measurements are used, and how their error prop- 
agate to the angular diameter distance. In section 3 we 
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give a detailed analysis of the error budget, in section 4 
we discuss the constraints on cosmological parameters. Fi- 
nally, section 5 summarizes our conclusions. 

2. DETERMINATION OF ANGULAR DIAMETER DISTANCE 
USING CLUSTERS 

Distance determinations using the SZ effect and X-ray 
emission from the intra-cluster medium (SZ/X method 
hereafter) are based on the fact that these processes de- 
pend on different combinations of physical parameters of 
the clusters. The SZ effect (Sunyaev and Zel'dovich 1980; 
for recent reviews cf. Birkinshaw 1999, and Rephaeli 1995) 
is a result of inverse Compton scattering of CMB photons 
off hot electrons in the intra-cluster (IC) gas. The number 
of photons is conserved, but, on average, the photons gain 
energy, and thus generate a decrement in the Rayleigh- 
Jeans part of the spectrum and an increment in the Wien 
region. The amplitude of the SZ effect does not depend 
on the redshift of the cluster. We will discuss static ther- 
mal and kinematic thermal SZ effects in this paper, where 
the "static" effect is present in all clusters, and the "kine- 
matic" effect is only present for those clusters with a non- 
zero line-of-sight (LOS) peculiar velocity relative to the 
Hubble flow. A typical Rayleigh-Jeans decrement of the 
static SZ effect is about about 50 times that of the kine- 
matic SZ effect. The static SZ effect is proportional to the 
LOS pressure integral of the IC gas 

(1) 



AT oc g{u) J dz n e (r) T e (r), 



where g(v) is the frequency dependence of the effect (g — > 
—2 in the non-relativistic Rayleigh-Jeans limit), and n e (r) 
and T e (r) arc the electron density and temperature as 
functions of position within the cluster. The central X-ray 
surface brightness of a cluster is emission- weighted line-of- 
sight average of n e (z) 2 T e {z): 

S * a n 1 m I dz n ^ z f T ^ z ) A ( T e> Z "b), (2) 
(1 + z) J 

where A(T e ,Z a b) is the cooling function integrated over 
the de-redshifted energy band of observations and Z a b is 
the metal abundance of the gas. 

On any given LOS, Equations l_and 2, and the emis- 
sion weighted X-ray temperature, T e (from spectroscopy), 
provide three independent integral constraints for the two 
functions, n e (rLOs) and T e {vLOs), although these func- 
tions cannot be determined uniquely. If, instead, one as- 
sumes parameterized functional forms for n e (r) and T e (r), 
these three equations can be used to constrain the control- 
ling parameters. One important parameter that can be 
found from fitting the data is a characteristic LOS phys- 
ical size of the system, R^hif' an d * ne angular diameter 
distance of the cluster can then be determined if this size 
is compared with the corresponding angular size, O^f. 
However, we can measure only the apparent characteristic 
angular size, 9f h K ^ of the system as it appears projected 
into the plane of the sky (with corresponding physical size 
R fhar)- Tnc main difficulty in the SZ/X method is the 
de-projection of the cluster from its 2-D images, that is 
the problem of finding the value of B^har lrom f ne mea- 
sured value OfhJr- If this can be done, one can determine 
the angular diameter distance as 
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If one assumes spherical symmetry, 0f^f = O^f, and 
a de-projection is possible without further assumptions 
about the functional forms of n e (r) and T e (r) (Silk & 
White 1978). This method, however, needs high sig- 
nal/noise and high angular resolution SZ and X-ray images 
of the cluster, and thus the analytic de-projection has not 
been used so far. Most commonly, the distribution of gas 
in clusters is described by the spherical isothermal beta 
model (Cavaliere & Fusco-Femiano 1976). This model 
seems to give a good fit to X-ray data (Sarazin 1988), and 
thus it is usually assumed that the IC gas follows a spheri- 
cal isothermal (5 model, where the IC gas is isothermal and 
in hydrostatic equilibrium in the gravitational potential of 
the cluster. The electron concentration in the cluster at- 
mosphere is then n e (r) oc (l+(r/r c ) 2 )~ 3/3//2 , where r c is the 
core radius (a suitable characteristic size of the cluster). 
The shape parameter j3 describes how the kinetic energy is 
distributed between the galaxies and IC gas. Typically r c 
= 0.1 Ii^qq Mpc, and (3 w 2/3 (Sarazin 1988). In practice, 
the IC gas temperature is determined from X-ray spec- 
troscopy, the core radius projected in the plane of the sky, 
0^ IT , and the shape parameter (3 FIT are determined from 
fitting the /? model to a high resolution X-ray image. Fi- 
nally equations 1 and 2 are used to determine the angular 
diameter distance using Equation 3 
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which is crucially dependent on the assumption of spheri- 
cal symmetry: rf IT — r^ os . 

If we adopt the beta model representation, we can write 
the angular diameter distance as a function of measurable 
quantities for a cluster as 

1 (AT ) 2 K{T e ,Z ab ) F(f3) 



D A oc 



(5) 



(3) 



(1 + z y S X o(N H ) t 2 9 2 Mfl c ' 

where ATo and Sxo are the central SZ effect and X-ray 
surface brightness (un-absorbed), Nu, is the absorbing 
column density, and F(/3) is a known function of (3 (or 
other parameters if a different description of the structure 
of the cluster is used). Some of the measured (or esti- 
mated) quantities on the RHS of this equation are highly 
correlated, so a careful error analysis is needed: the er- 
rors on individual quantities cannot be added in quadra- 
ture. No bootstrapping is required in using Eq. 5 (it is 
a direct method, with no need for using the usual cosmic 
distance ladder) as was realized by Sunyaev & Zel'dovich 
(1980), but the strong assumptions of spherical symmetry 
and isothermality are only crude approximations to the 
real situation. 

It has been pointed out recently that relativistic effects 
become important for high temperature clusters (Rephaeli 
1995). Relativistic corrections to the inverse Compton 
scattering have been calculated and used to interpret ob- 
servations in terms of the Hubble constant (Rephaeli 1995; 
Birkinshaw 1999; Birkinshaw & Hughes 1998; Holzapfel 
et al. 1997, and references therein). Relativistic cor- 
rections to the thermal bremsstrahlung cooling function 
are also needed (Hughes & Birkinshaw 1998; Rephaeli & 
Yankovitch 1997). These corrections are in the order of 
3% in Da for high temperature (T e > 10 keV) cluster. 
In this paper we assume that relativistic corrections have 
been made and we will not discuss their effect in detail. 



4 



3. ERROR BUDGET IN THE DETERMINATION OF THE 
ANGULAR DIAMETER DISTANCE 

Errors in the determination of the angular diameter 
distance may be cast into two major categories: errors 
from measurements, and errors from theoretical modeling. 
Measurement errors can be statistical or systematic, while 
errors from theoretical modeling are systematic by nature. 
The systematic errors can be further classified as "ran- 
dom" or "non-random" depending on whether they aver- 
age out or not for a statistical sample of clusters. In what 
follows we therefore use the expressions "random modeling 
errors" when discussing errors which introduce only scat- 
ter in the distance determination when averaged over an 
unbiased sample of clusters, and "non-random modeling 
errors" when discussing errors which introduce a bias in 
the distance determinations even for an unbiased sample 
of clusters. 

In this Section we summarize the statistical and system- 
atic errors associated with the SZ and X-ray observations, 
following the discussions in reports of the latest interfero- 
metric observations (from BIMA, Reese et al. 2000; and 
from the Rylc Telescope [RT], Grainge et al. 1999), and in 
review articles (Birkinshaw 1999). We also use some addi- 
tional references to compile a detailed list of error sources. 

3.1. Errors from measurements 

Statistical errors from measurements in the angular di- 
ameter distance are dominated by counting statistics in 
X-ray images and spectra, and by Gaussian measurement 
uncertainties in the SZ measurements. These statistical 
errors propagate to errors in 9 C and (3 through fitting for 
the spatial distribution, and to errors in T e , A e , Z a b, and 
Njj through fitting the X-ray spectrum. A simultaneous 
fit of a (3 model on interferometric SZ and X-ray images 
determines 6 C , (3, Sxo and AT (e.g. Reese et al. 2000), 
but causes the errors in these parameters to be strongly 
correlated. 

The most important error sources in a determination of 
the angular diameter distance are the error in Sxo and 
the error in T e (Eq. 5). The electron temperature enters 
the estimate of Da roughly as XjT 2 , since for most X-ray 
observations to date the X-ray emissivity, A e (T e , Z a b) is a 
slowly- varying function of T e . Other uncertainties, from 
the dependence of A e on the abundance of metals in the 
ICM and the absorbing column on the line of sight, are less 
important, and we can also neglect the contribution to the 
error in the angular diameter distance from the redshift. 

The error in the spatial fit is dominated by the uncer- 
tainty in the central SZ decrement, ATq, and central X-ray 
surface brightness, Sxo- Both errors are about 8-10% in 
the data obtained using BIMA and ROSAT. The total 
uncertainty in a Da estimate from spatial fitting to X- 
ray imaging and SZ interferometric measurements (with 
BIMA or the RT) is about 14-18%. We can expect a dra- 
matic improvement in the accuracy of Sxo i n the near 
future because of the larger collecting area and higher an- 
gular resolution of XMM and Chandra, and because the 
improved imaging will also allow a better choice of models 
for the ICM. A substantial improvement in the central SZ 
effect is also likely, as the first generation of dedicated SZ 
interferometers becomes available. The overall statistical 
error in these parameters should drop to about 4-5%. 



The statistical errors in measurements of electron tem- 
perature based on ASCA and ROSAT observations are 
up to about 15-20% (Reese et al. 2000; Grainge et al. 
1999). With today's technology (Chandra and XMM), a 
5% statistical error in the emission-weighted average elec- 
tron temperature of a well-defined region of a high redshift 
cluster should be straightforward (David et al. 2001; Ar- 
naud et al. 2001; Peterson et al. 2001). 

The principal identified systematic errors in the mea- 
surements arise from the absolute calibration of the radio 
and X-ray observations. Errors in the effective area of the 
ROSAT PSPC and HRI introduce an error of about 10% 
in Da- This error is greatly reduced (to 1-2%) for the in- 
struments on Chandra and XMM-Ncwton. The absolute 
calibration error of radio interferometers is good to 4-5%, 
but could be improved to about 1% through extensive ob- 
servations of point sources, tied to the planetary flux den- 
sity/brightness temperature scale (Birkinshaw 1999). 

Ground-based interferometric observations are also sub- 
ject to systematic errors due to the removal of background 
(and cluster) radio sources, which may be imperfect if the 
sources are variable or have significant angular extent. Im- 
perfect calibration of the phase and amplitude of the de- 
tector system may also cause errors near brighter contami- 
nating sources. Fortunately interferometers are insensitive 
to large scale gradients in emission from the ground and 
atmosphere, and so this source of systematic offset signals 
from single-dish data is largely removed in interferometric 
work (e.g., Carlstrom et al. 2001). 

3.2. Errors from modeling 

The modeling of the structure of the ICM is an impor- 
tant part of determining the distance to a cluster, since it 
is the model that makes it possible to connect the LOS size 
of a cluster to its apparent angular size projected on the 
plane of the sky. Following most treatments in the litera- 
ture, we use isothermal spherical (3 models to describe the 
IC gas distribution, and in this Section discuss the errors 
introduced by deviations from this model. As before, we 
can identify random and non-random systematic modeling 
errors. 

Significant random modeling errors arise from the as- 
phericity of the intra-cluster gas, the peculiar velocity of 
the cluster, and primordial CMB fluctuations. Serious 
non-random modeling errors can be expected from non- 
isothermality, cooling flows, clumping, merging, and finite 
extent of the cluster. Other issues that may be impor- 
tant are radio point sources and gravitational lensing. We 
briefly discuss random modeling errors from resolved radio 
halos in clusters because of their contribution to the mea- 
sured SZ signal (through synchrotron emission) and to the 
X-ray signal (through inverse Compton emission), the ef- 
fect of diffuse free-free emission from cool gas (for example 
in spiral galaxies), and finite optical depth effects. 
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3.2.1. Random modeling errors 

Asphericity is one of the most important source of sys- 
tematic errors in the determination of the distance to 
clusters when using a spherically symmetric isothermal (5 
model. If the cluster is not spherical, the assumption that 
r c IT = f"c 1S invalid. We estimate the extent of this 
problem by approximating the true structures of clusters 
as oblate or prolate ellipsoids (or, more generally, as tri- 
axial ellipsoids), while retaining the assumption that the 
distribution is described by an isothermal /3 model with 
constant [3. If the symmetry axis is in the LOS, we would 
assume r^ IT — r^ os , and thus overestimate (underesti- 
mate) Da for prolate (oblate) clusters. Birkinshaw (1999) 
finds that prolate or oblate clusters with symmetry axis 
aligned in the LOS, if assumed spherical, will have a frac- 
tional error of 



where a(= r^ os ) is the core radius in the LOS and 
b(= r^ IT ) is the core radius in the plane of the sky (a > b 
for a prolate distribution). Clusters often show elliptic- 
ity at the level, in projection, of a/b = 1.25 (Mohr et 
al. 1995). If the axis of symmetry is not in the LOS, 
the error is smaller, therefore the fractional error in Da 
should be < 25%. Hughes & Birkinshaw (1998)'s analysis 
of CL0016+16 showed that oblate or prolate distributions 
may cause less than 8% error in Da if the structure of 
the ICM is analyzed as spherical. Grainge et al. (1999) 's 
work implies about a 14% error in Da for Abell 1413. 

N-body simulations can be used to understand the de- 
tails of physics of the cluster geometry, and quantify the 
deviations from the spherical distribution. Inagaki et 
al. (1995)'s numerical simulations show that aspheric- 
ity causes an error of up to 15% in Da- Cluster merging 
simulations of Roettiger et al. (1997) show that small 
off axis merging in general causes prolate distributions, 
while oblate distributions may be caused by large off axis 
merging, causing typically less than 20% error in Da- Ob- 
servationally, Basilakos et al. (2000) found that prolate 
spheroid models fit the APM cluster data better than 
oblate spheroids. Sulkanen (1999) studied a statistical 
sample of clusters assuming a triaxial beta model density 
distribution. Sulkanen's results indicate that the distance 
scale obtained assuming a spherical distribution is within 
14% of its true value (at 99.7% confidence) based on a 
sample of 25 clusters with triaxial axes consistent with 
observations (cf. also Puy et al. 2000). 

In general, oblate or prolate clusters will have their axis 
randomly distributed in the sky. Zaroubi et al. (1998) 
studied general de-projections assuming that clusters have 
axially-symmetric density distributions with an arbitrary 
orientation of the symmetry axis. They found that using 
SZ and X-ray images one can determine only a combi- 
nation of distance and inclination angle. Another image 
(weak lensing for example) is necessary to decouple these 
two parameters and determine the distance separately. If 
only X-ray and SZ data are available, the error that as- 
phericity introduces into the Da determination for any 
single cluster cannot be reduced below about 15%. How- 
ever, the combination of X-ray, SZ and weak lensing data 
together with modeling of the equilibrium of the hot gas, 
should allow us to reduce the error to around 5%. The 



planned weak lensing surveys should make the required 
lensing data available in the near future. Alternatively, 
galaxy velocity distributions from optical observations can 
be used with numerical simulations to determine physical 
parameters of individual clusters, as was done by Gomez et 
al. (2000) and Roettiger et al. (1997) (see the discussion 
of merging, below). 

Peculiar velocities of clusters introduce enhanced (ap- 
proaching) or decreased (receding cluster) SZ measure- 
ments because of the kinematic SZ effect (Sunyaev & 
Zel'dovich 1980; for a derivation cf. Birkinshaw 1999). 
In CDM models, cluster peculiar velocities are 400- 
500 kms" 1 (Colberg et al. 2000; Ueda et al. 1993). This 
would introduce only a few percent error in Da- How- 
ever, some observations suggest larger peculiar velocities, 
1000 kms -1 (Bahcall & Soneira 1983; Lauer & Postman 
1994). If these large velocities are real, the kinematic SZ 
effect may cause an error of up to 25% in Da- Fortunately 
one can separate the static and kinematic SZ effects based 
on their different frequency dependence or, equivalently, 
the peculiar velocities can be determined by measuring 
the cross-over frequency of the total SZ effect (e.g., Mol- 
nar & Birkinshaw 1999). Peculiar velocities also introduce 
a small bias in the redshift determination of the cluster. 
If clusters are selected based on their SZ signal, this effect 
would cause a biased sample of clusters, and a systematic 
overestimate of Da would result. 

Primordial CMB fluctuations introduce systematic ef- 
fects in the distance determinations with their positive or 
negative contributions to the microwave decrement mis- 
interpreted as SZ signal. The amplitude of systematic 
errors introduced by CMB fluctuations is a strong func- 
tion of the observation strategy. The CMB fluctuations 
are reduced at arcminute angular scales compared to their 
degree-scale amplitudes, but still introduce a scatter of 
about 10% in the distance determinations (Cen 1998). 
At smaller scales (high £s) the power in CMB fluctuations 
becomes negligible. A further reduction in the level of 
this error can be achieved using spectral separation of the 
thermal SZ effect from the primordial fluctuations (and 
the kinematic SZ effect). Gravitational lensing transfers 
power from large scale primordial fluctuations to small 
scale fluctuations (Metcalf & Silk 1998; Seljak 1996). As 
a result, if the CMB is not separated from the SZ effect, 
it would give larger, but a symmetric scatter in Da, or a 
fractional error of about 8% (Cen 1998). Note, however, 
that the kinematic SZ effect and primordial CMB fluctua- 
tions have the same frequency dependence, thus they can 
not be separated from each other based on their frequency 
signature. Fortunately both effects can be separated out 
from the static SZ effect simultaneously, based on their 
different frequency dependence. 

The errors in the cosmological parameters caused by 
random modeling errors (asphericity, cluster peculiar ve- 
locities, and primordial fluctuations, etc..) can be reduced 
using a properly-selected sample of clusters. This sample 
must avoid selection biases that themselves introduce sys- 
tematic errors. As was emphasized by (Birkinshaw et al. 
1991), clusters should not be selected based on their SZ 
or X-ray central brightnesses, since such a selection would 
produce a sample containing an excess of clusters with pro- 
late geometry, high positive peculiar velocity in the LOS, 
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and contamination from negative CMB fluctuations (if the 
Raylcigh- Jeans frequency band is used), and result in a 
biased Da- X-ray selected clusters with a flux limit well 
above the detection limit might be used as in Mason et 
al. (2001) and Jones et al. (2001). Alternatively, if 
there are multi-frequency measurements, the best solution 
would be to separate the kinematic SZ effect and the pri- 
mordial fluctuations from the static SZ effect, and use the 
total static SZ effect flux density as a selection criterion. 

3.2.2. Non-random modeling errors 

Non-isothermality is one of the major sources of system- 
atic error in the determination of Da- 

The isothermal assumption for the intra-cluster gas is 
clearly an approximation. Even if the central region of 
cluster is virialized and isothermal, the outer regions will 
be subject to shocks from merging and gas in-falling from 
filaments. Merging with massive clusters will change the 
temperature relative to the single-component cluster viri- 
alized value, even in the core region. Observations show 
that the central regions are nearly isothermal, but thermal 
substructures have also been found in clusters (Sarazin 
1988). The effect of temperature variations on the dis- 
tance estimate depends on the instrument and observing 
technique used. SZ measurements are insensitive to tem- 
perature variations in the cluster if the projected pres- 
sure profile is unchanged. Non-spatially resolved X-ray 
measurements determine emission weighted temperatures, 
and so are sensitive to the thermal structure of the central 
region of clusters. Thus the temperature deduced from 
X-ray measurements is well suited to comparison with SZ 
measurements by radio interferometers, which also most 
sensitive to the central region and have less response to 
the outer parts of the cluster (for example, assuming 2a 
detection, BIMA is sensitive out to about 3r c , and about 
85% of the observable X-ray flux of a typical cluster, using 
ROSAT PSPC, is within 3r c ). 

Birkinshaw & Hughes (1994) and Holzapfel et al. 
(1997) analyzed Abell 2218 and Abell 2163 respectively, 
assuming an isothermal /? model and a model with falling 
temperature with radius. Their results show that Da may 
be overestimated by 20-30% if non-isothermal distribu- 
tions are assumed to be isothermal. Numerical simula- 
tions show similar errors due to non-isothermality: thus 
Inagaki et al. (1995)'s simulations lead us to conclude 
that an overestimate of 25% in Da, mostly because of the 
overestimated SZ amplitude, may result from assuming 
isothermality when the temperature is lower in the outer 
regions. Simulations of merging clusters by Roettiger et 
al. (1997) obtained a similar result, that an overestimate 
of 10-30% may result in Da from non-isothermality, with 
the range depending on the projection geometry. Note, 
however, that these results are based on single-dish mea- 
surements: for interferometric observations these effects 
are usually smaller. 

It is not easy to correct for non-isothermality. The SZ 
effect is proportional to < n e T e > R^haf' anc ^ f° r an ac ~ 
curate calculation of the SZ effect, it is often necessary to 
have good information about the temperature out to R V i r 
(about 10 r c ) or more. It is difficult to carry out spatially 
resolved spectroscopy at such large radii because of the 
low X-ray surface brightness of the outer regions of clusters 



(Sx oc rig). At present, the best evidence from BcppoSAX 
data extends T(r) measurements to (0.5 — 0.75)R V i r (Ir- 
win & Bregman 2000; De Grandi & Molendi 1999). The 
increased sensitivity available with Chandra and XMM 
should enable us to use spatially resolved spectroscopy to 
determine the temperature profiles of clusters which both 
resolve the inner cooling flow region and collect enough 
photons to measure useful temperatures in the clusters' 
outer regions (Schmidt, Allen & Fabian 2001; Tamura et 
al. 2001). 

Unfortunately the hydrostatic equilibrium assumption 
of the (3 model breaks down in the central parts of clus- 
ters with high enough density for cooling to be impor- 
tant. In these cases the ICM will radiate via thermal 
brcmsstrahlung and line emission (with the balance de- 
pending on the temperature), and develop a pressure gra- 
dient and a sub-sonic inflow, a so-called cooling flow (see 
review of Fabian 1994). The increased central density at 
the core of the cluster leads to an increased level of X- 
ray emission, which is often used as an indicator of the 
presence of a cooling flow region. 

Phenomenological models have been developed and nu- 
merical simulations have been performed to study cooling 
flows (White & Sarazin 1987; Rizza et al. 2000; Fabian 
1994, and references therein). Majumdar and Nath (2000) 
estimated the effect of cooling flows on the determination 
of the Hubble constant. They found that only at the very 
center of the cooling flow (within the sonic radius) will 
the X-ray luminosity drop because of the decreased tem- 
perature of the ICM. They show, further, that there is 
a gradual increase in pressure towards the center of the 
cooling flow. Outside this region, the approximation of 
hydrostatic equilibrium profile remains good, while within 
there will be an SZ effect excess. Thus Majumdar & Nath 
(2000) 's results indicate that one should exclude 80% of 
the cooling flow region to reduce the error in Da to below 
10%. However, this may be an upper limit on the error in 
Da, since the calculation ignored the effect of the cooling 
flow on the structure fitting. Fig. 3 and 4 in Reese et al 
(2000) suggest that a compensatory error occurs here, and 
this lowers the error in Da- 

A further effect from cooling flows, pointed out by 
Schlickciscr (1991), is that the build-up of cold gas at the 
center of cooling flows might lead to significant free-free 
emission in the radio band, which would reduce the SZ 
signal. This works in the opposite sense to the pressure 
effect, but provides another reason for excluding the cen- 
tral parts of SZ images of cooling flow clusters from the 
SZ/X distance scale analysis. 

Excluding cooling flow clusters completely from SZ/X 
distance scale studies would be the most complete solu- 
tion to the problem that they pose. However, most clus- 
ters close to hydrostatic equilibrium (where the underlying 
physics required to model the ICM is most straightfor- 
ward) possess cooling flows. Modeling cooling flows, and 
excluding their centers, is therefore necessary to build up 
large samples for Hubble constant work. 

High-redshift cooling flow clusters can be recognized by 
using their X-ray spectra even if we cannot resolve the 
central region, since cooling flow clusters have emission 
weighted metallicity 1.8 times higher than non-cooling flow 
clusters (Allen & Fabian 1998). Mohr, Mathiesen & 
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Evrard (1999) analyzed nearby cluster of galaxies using 
ROSAT PSPC data, modeling cooling flow regions where 
necessary, and found that fitting an isothermal (3 model to 
cooling flow clusters will underestimate both r c and /3, and 
therefore produce a poor fit even outside the cooling flow. 
They identified cooling flow clusters based on two criteria: 
1, non-random residuals consistent with a central emission 
excess; and 2, relaxed cluster with no asphcricity or sub- 
structure. They concluded that a double (3 model fit gives 
an unbiased estimate of Da- Cooling flow contamination 
in the determination of the average Hubble constant using 
the SZ/X ray method can be recognized by searching for 
a dependence of individual Hubble constants (determined 
for each cluster) on IC gas metallicity. 

Clumping in the IC gas is potentially one of the most im- 
portant systematic effects. It is well known that radio ha- 
los, hot bubbles from supernova eruptions, cold condensed 
gas in galaxies, etc. constitute a level of dumpiness in 
the ICM. The important question is whether the effect of 
clumping is important in the SZ/X method of measuring 
D A - 

As a first approximation, clumping will enhance the X- 
ray surface brightness since Sx cx n 2 , and does not change 
ATsz, since AT$z is proportional to the line-of-sight av- 
eraged pressure and clumps should be in pressure equilib- 
rium with their surroundings if they are to be long-lived 
(assuming no substantial magnetic fields exist) . As a con- 
sequence, Da is underestimated. However, the emission- 
weighted temperature of the cluster measured by X-ray 
spectroscopy will also decrease. This reduces the system- 
atic effects of clumping to only a few percent error in Da- 
Birkinshaw, Hughes & Arnaud (1991) studied the effects 
of isobaric clumping of the intra-cluster gas in Abell 665. 
The fractional error in Da from clumping is 

SDa <n e > 2 < A(T e ) > _ 

D A < n 2 A(T e ) > ' U 

where the angle brackets imply averages over regions larger 
than the scale of the clumping. For isothermal clumps 
A e factors out. Then since < n\ > / < n e > 2 is al- 
ways greater or equal to 1, clumping will always cause 
an underestimate of Da- Birkinshaw et al. (1991) found 
< n\ > 2 / < n e > 2 < 3 or so. Holzapfel et al. (1997) studied 
the effects of isobaric clumping on Abell 2163. Assuming 
cold clumps they estimate that error from clumping could 
lead to an overestimate of about 10% in Da- 

Inagaki et al. (1995) find, from their fl m — 1 numeri- 
cal simulations, that Da will be underestimated by about 
15% due to clumping. However there is evidence that lower 
matter density models produce less clumpy structures, so 
that if our Universe has Q, m < 1, Inagaki et al. will have 
overestimated the effect of clumping. 

Self-consistent modeling of small scale clumping is diffi- 
cult because many physical processes contribute to its cre- 
ation and destruction. Gunn & Thomas (1996) used phc- 
nomcnological multi-phase models to study X-ray emission 
from clumpy IC gas. Their isobaric model implies a frac- 
tional error in the X-ray central surface brightness 
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Nagai et al. (2000) discussed biases in the Hubble 
constant determination from a multi-phase, spherically- 
symmetric, intra-cluster medium with isobaric clumping 
of variance 

° 2 = (l + (r/r c ) 2 )^' (9) 

where a c and e are free parameters describing the strength 
and radial dependence of the variance. They assumed a 
log-normal distribution for the gas density phase distribu- 
tion, which was motivated by simplicity and the effect of 
non-linear gravitational growth of initially Gaussian den- 
sity fluctuations (Cole, Fisher and Weinberg 1994). They 
assumed that the multi-phase model has the same emis- 
sion weighted temperature and the X-ray emission profile 
as a fiducial single phase model. Based on their results, 
the error in the distance from an incorrect assumption of 
a single-phase ICM is 



^ = 2cx P [iI^^ 2 ] 
Da 4 



(10) 



where a is the emissivity exponent. This would cause 
about a 10% error in Da- 



where a is again the power-law exponent in the emissiv- 
ity function. We can conclude that clumping may cause a 

5 — 20% error in Da- 

Spatially-resolved X-ray spectroscopy, as can be per- 
formed using Chandra and XMM, will help to estimate the 
dumpiness of the IC gas. Emission lines between 0.5 and 
1.5 keV originating in cool regions, such as the Fc L-shell 
lines, H- and He-like lines from N, O, Ne, Mg, will provide 
a strong test on multi-phase models. Unfortunately this 
method will not constrain all types of dumpiness. Hughes 

6 Birkinshaw (1998), and Mason (1999) suggest that, 
since SZ/X Hubble constant determinations are not very 
different from the results of other methods, clumping in- 
troduces less than a 10% error. However, it is possible that 
clumping itself introduces more bias, but that this bias is 
balanced by other effects with opposing biases. 

Roettiger et al. (1997) studied systematic errors in the 
Hubble constant from cluster merging. Merging can lead 
to the formation of shocks that compress and heat the 
IC gas. After a merger, the gas will settle into the pro- 
late or oblate potential well defined by the dark matter 
distribution. As already discussed, ellipticity introduces 
scatter rather than bias in Da- However, shock- heated 
IC gas will have different SZ and X-ray properties than 
gas in hydrostatic equilibrium, and its presence will intro- 
duce bias. Based on simulated X-ray images and tempera- 
ture maps, Roettinger et al. find that isothermal assump- 
tions in merging clusters systematically overestimate Da 
by about 15%. Numerical simulations show that density 
enhancements due to merging of sub-clusters can result in 
higher X-ray surface brightness, change mass estimates, 
and might, if not recognized, cause an about 20% error in 
distance determination (Mohr et al. 1999). 

Roettiger et al. (1997) suggest that clusters at the 
early stages of merging should be excluded from distance 
determinations. Dynamically active clusters may be rec- 
ognized from their galaxy velocity distributions. Clusters 
with dynamical activity will have a large (3 discrepancy 
Wfit 7^ Pspect, see details in Sarazin 1988). Anisotropy in 
galaxy velocity distribution also signals dynamical activ- 
ity. More relaxed systems, like merging clusters after they 
reached quasi-equilibrium may be modeled and included in 
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the distance determination. Simulations help to analyze 
individual merging clusters (Roettiger et al. 1995; Gomez 
et al. 2000): adjustments to the initial parameters of 
merging clusters are made until the resulting merged clus- 
ter has the observed SZ and X-ray appearance and galaxy 
velocity distribution. 

The (3 model gives divergent masses if not truncated at 
some finite radius, which is usually taken to be about 10 
core radii. This fact is a sign that at large radii, the beta 
model can not be correct. When calculating the distance 
to a cluster using the isothermal (3 model, we assume an 
infinite extent, and therefore the finite extent of clusters 
introduces a systematic bias in Da- The SZ effect is more 
sensitive to outer regions than X-ray bremsstrahlung, since 
it is proportional to n e unlike X-ray bremsstrahlung, which 
is prop to n\. In theory, the SZ effect measurements are 
more suitable for studying the outer regions of clusters, 
and also are subject of more error caused by the finite 
extent of clusters. 

Inagaki et al. (1995) discussed the effect of finite clus- 
ter sizes. They find that the ratio between finite truncated 
and full beta models is 

r „ r, 5,(3^/2-1/2,1/2) -I 2 

dL> A = 5(3/3/2-1/2,1/2) , . 

D A 1 - (3/3-1/2,1/2) ' ^ I 

L 1 3(3/3-1/2,1/2) J 

where q = 1/(1 + p 2 ), p = i? cut /r core , B and B q are 
the beta and incomplete beta functions, and assuming no 
change in 9^ IT . Inagaki et al. (1995)'s results suggest 
that this effect can cause an about 10-20% underestima- 
tion of Da for typical parameter values (cf. also Puy et 
al. 2000). By contrast, Birkinshaw & Hughes (1994) 
and Holzapfel et al. (1997), analyzing Abell 2218 and 
Abcll 2163 respectively, find that finite cluster extent con- 
tributes only about 6% and 2% overestimate error in Da- 
Observationally, Molnar (2000) attempted to find the 
outer cut off of Abell 3571 using RXTE scans across the 
cluster. Abell 3571 seems to extend well beyond its virial 
radius, however, this conclusion is not strong due to poor 
statistics in the data. 

Clearly, more observations are needed to find out the 
extent to which the (3 model is correct, and how the tran- 
sition happens from the cluster to the surrounding regions 
(dominated by filaments according to numerical simula- 
tions). A falling temperature profile with radius is a sign 
of deviations from the (3 model, and may provide the best 
tracer of additional structure of this type. We expect er- 
rors from this effect to be reduced greatly when the new 
spectral information from Chandra and XMM is available. 

One of the most important contaminants of the SZ ef- 
fect is emission from radio point sources. Such emission, 
especially at the center of the cluster, will decrease the fit- 
ted amplitude of the SZ decrement in the Rayleigh-Jeans 
frequency region, and cause an underestimate in the dis- 
tance of the cluster. Using the sensitivity of the instru- 
ment, one can estimate the maximum flux density of un- 
resolved point sources that can contribute to this error: for 
BIMA and the RT, the systematic errors are in the 10-15% 
range. The level of radio source confusion varies with the 
observing technique used. Multi-baseline interferometric 
observations have the most favorable confusion level. A 
deep survey, carried out at different frequencies or inter- 
ferometer baselines, would help to find point sources up 



to a limit when the error due to unresolved point sources 
would be negligible relative to other errors. Interferomet- 
ric measurements also have the advantage over single-dish 
observations of being able to monitor the brightnesses of 
(potentially-variable) point sources while simultaneously 
measuring the SZ effect. 

As Loeb & Refregier (1997) pointed out, systematic ef- 
fects in the determination of the baseline for the SZ effect 
arise from gravitational tensing by the cluster gravitational 
potential. The brightnesses of point sources are enhanced 
by lensing, which brings them above the detection thresh- 
old, and thus they are removed from the field. This over- 
removal of point sources from the background lowers the 
background flux relative to a control field, and thus leads 
to an overestimate of the SZ signal, and an underestimate 
of about 10% of Da- However, this effect is important 
only at frequencies less than about 30 GHz, and its pres- 
ence can be tested (and corrected for) using the model 
cluster mass distribution which can be obtained by con- 
ventional analyses of the X-ray data for each cluster (cf. 
also Blain 1998). 

Other, probably small, effects contributing to the er- 
ror budget, which should be checked and treated indi- 
vidually when necessary are: contributions in the radio 
band from synchrotron emission from resolved halo (and 
other) sources in clusters, free-free emission from cool gas 
in spirals, free- free emission from radio halos in the X- 
ray band (Birkinshaw 1979), and any contribution from 
the non-thermal SZ effect. As was pointed out by Mol- 
nar & Birkinshaw (1999), the Kompaneets equation and 
its relativistic extensions are equivalent to a single scat- 
tering approximation, thus a small effect will arise from 
ignoring finite optical depth (multiple scattering). Any 
non-thermal population of IC electrons would produce a 
non-thermal SZ effect. There are some theoretical argu- 
ments for such a population (Petrosian 2001; Blasi 2000; 
Colafrancesco 1999; Sarazin 1999), and some observa- 
tions have detected the expected excess hard X-ray emis- 
sion (e.g., in the Coma cluster), but such excess X-ray 
emission is rare and hence the non-thermal electron pop- 
ulation is weak (Fusco-Femiano et al. 2001; Maloney & 
Bland-Hawthorn 2001). Clusters with complex morphol- 
ogy should not be used for distance determinations be- 
cause of the difficulty in their modeling. If substructure in 
high redshift clusters (as in RX J1347-1145, cf. Komatsu 
et al. 2001) is common, one should check such individual 
clusters for complex morphology using a high resolution 
instrument. 

The cumulative systematic effect resulting from all these 
sources of error is best addressed via numerical simula- 
tions. Yoshikawa et al. (1998) used numerical simulations 
to estimate the systematic error from fitting (3 models to 
clusters. They used a spatially flat fiducial model with 
Vl m = 0.3, f^A = 0.7, and h = 0.7 for their simulations. 
They found that the isothermal (3 model describes clusters 
well both in SZ and X-ray imaging. They found that, even 
though the fitted (3 model parameters were different when 
fitted to SZ or X-ray images (due to non-isothermality, 
non-asphericity, and dumpiness), the systematic errors 
in the Hubble constant (and so in distances) are negli- 
gible at low rcdshifts. At high redshift (z ~ 1), about 
a 20% overestimate of Da might occur, mainly due to 
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non-isothermality and asphericity. Inagaki et al. (1995) 
find that the cumulative effect of non-isothermality and 
asphericity leads to about a 10 — 20% overestimate of Da, 
since the overestimate due to non-isothermality is larger 
than the underestimate due to dumpiness (though in their 
simulations fl m = 1). 

3.3. Summary of the error budget 

The new generation of satellites will allow us to deter- 
mine more precisely the temperature profile and spatial 
structure of the intra-clustcr medium, and thus minimize 
the systematic errors in Da that result from temperature 
variations. When necessary, the spatial resolution of these 
new instruments will allow us to model the gas in each indi- 
vidual cluster beyond the spherical isothermal beta model. 
Numerical modeling of individual clusters will help us to 
derive their physical parameters. The improved spatial 
resolution will also allow us to study dumpiness, which 
is another important source of uncertainty. In general, 
those effects which have different spectral signature from 
the static SZ effect, should be separated from the SZ effect 
using multi- frequency observations. 

Based on our evaluation of the error budget, we esti- 
mate that the random error in Da achievable in the near 
future using known techniques, and assuming that lensing 
measurements will be used to eliminate errors from clus- 
ter asphericity, might be as low as 7%. This estimate is 
made up by quadrature combination of the errors, with 
the dominant terms being from T e (5%) and spatial fitting 
(5%). A systematic error of 5% might also be obtained, 
although this would be difficult. 

Evolution effects are clearly important limitations in de- 
termining cosmological parameters. The advantage of the 
SZ/X-ray method is that, as long as hydrostatic equilib- 
rium and simple geometry hold, it should be reliable even 
if scaling laws (for example mass - temperature) evolve. 
However, a detailed analysis of evolution effects is out of 
the scope of our paper. We restrict ourselves to simply 
studying the effect of an additional systematic error in the 
Da determination with a linear gradient with redshift. 

4. CONSTRAINTS ON COSMOLOGICAL PARAMETERS 

As discussed earlier, in the near future SZ surveys will 
discover hundreds of high and low redshift clusters, and we 
can expect Chandra and XMM observations of hundreds 
of clusters. However, accurate angular diameter measure- 
ments require long SZ and X-ray integrations and therefore 
we do not expect all discovered clusters to have accurate 
SZ/X distance measurements. If the SZ/X method is to 
be used to measure cosmological parameters, we should 
select a sample of clusters that minimize the systematic 
errors (Section 3) while not requiring excessive observing 
time. As pointed out in Hutcrcr & Turner (2000), the 
ideal distribution of clusters in redshift space would be a 
superposition of P delta functions in redshift, where P is 
the number of cosmological parameters to be determined. 
In practice, however, in any narrow redshift band we will 
not have enough clusters to achieve good statistics, and 
objects at z > 1 would need excessive integration times 
to reach high individual accuracies. A detailed study to 
select the optimum redshift distribution of clusters, taking 
into account the differing numbers of suitable clusters at 



different redshifts and the detailed characteristics of the 
possible instruments and observational strategies is be- 
yond the scope of our paper and more properly devolves 
on the groups proposing to construct such instruments. 
Here we do a simpler problem, by assuming that distances 
with similar accuracies are available for a set of clusters 
uniformly distributed in redshift space between z = 0.01 
and z = 1. This choice covers almost half of the red- 
shift interval most sensitive to the cosmological constant 
(z = 0.5 — 1.8), but excludes the more distant objects for 
which excessive integration times would be needed. 

We carried out simulations to estimate the constraints 
from the angular diameter distance-redshift function on 
the parameter space defined by ft m , Q\, and h, and also 
by O m , w, and h (assuming a spatially flat geometry 
Oa = 1 — ^m)- We simulated clusters using a spatially 
flat fiducial CDM model with Q m = 0.3, Cl\ = 0.7, and h 
= 0.65. We use the x 2 statistic to evaluate errors from our 
simulations, since we are dealing with large errors, 50% in 
cosmological parameters, and therefore the assumptions 
leading to the Fisher matrix formalism are not satisfied. 
We choose realizations which have best fit values close to 
those of the input fiducial model, and offset by the values 
of Ax 2 = 3.53, 8.02 and 14.2 appropriate for the number 
of fitted parameters to find the Gaussian 1,2, and 3 a er- 
ror surfaces (defined by probability levels of 68%, 95.4%, 
and 99.73%). 

As a first realization, we assumed a sample of 500 clus- 
ters uniformly distributed in redshift (so that 250 lie at 
z > 0.5) with 4% random error in Da, which might be 
achieved if systematic errors can be tightly controlled (we 
assume lensing measurements will be used to eliminate er- 
rors from cluster shape). We use these results to demon- 
strate the ultimate constraints on cosmological parameters 
that might be achievable from the SZ/X-ray method. In 
Figures 2a, b, and c, we show the results of fitting cosmo- 
logical parameters to the observationally-determined an- 
gular diameter distance function (1, 2, and 3a concentric 
ellipses, solid lines, assuming three independent parame- 
ters, corresponding to A% 2 of 3.53, 8.02 and 14.2 appro- 
priate for three fitted parameters). 2-dimensional (2D) 
projections of the 3D surface of 3a constraints are shown 
using dotted lines. In the Q m - Qa plane (Figure 2a) the er- 
ror ellipses are elongated roughly along the line 2 Q m — Qa 
= constant. This is similar to the constraint obtained 
from the SNe la experiment, which measures cosmological 
parameters via the predicted apparent magnitude distri- 
bution, and uses different redshift limits. Very low red- 
shift measurements of clusters lead to error ellipses which 
are elongated along the line fl m — 2fl\ = constant. As 
higher redshift objects are added, this elongated feature 
rotates counter-clockwise. From Figure 2c, which shows 
constraints in the h - Qa plane, we can conclude that 
with the assumed accuracy, the cosmological constant is 
not well constrained. However, h is well measured even 
without additional information. In general, the disper- 
sion of measurements (1, 2, 3a contours) is determined by 
random and not systematic errors. Systematic errors in- 
troduce only bias, shifting the mean of the measurements 
away from the expected value in the parameter space while 
preserving their dispersion. A ±3% systematic error in the 
Da will not affect the results in the Q m — plane, but 
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simply shift the error ellipses we obtained from random er- 
rors up or down along the h axis by 0.02 (3% of the fiducial 
value of h), since the amplitude of Da is set by the Hubble 
constant (Figures 2ab and c, 3a solid ellipses above and 
below the 1, 2, and 3a concentric ellipses corresponding 
to random errors, solid lines). With the assumed random 
and systematic errors, the determination of h becomes lim- 
ited by systematic errors. We also carried out simulations 
assuming a systematic error with a gradient in redshift 
growing from 0% at z — to ±3% at z = 1 in addition 
to the assumed 4% random error to study the effect of 
evolution. (Figures 2b and c, short dashed and dash dot- 
ted 3cr lines) . The error ellipses from random errors will be 
shifted in the Cl m - ft a plane (Figure 2a). From Figures 2b 
and c we can conclude that the Hubble constant will not 
be affected by this type of systematic errors, which fol- 
lows from the fact that the Hubble constant is constrained 
by the low redshift regime where the assumed systematic 
errors are small. f2 m and Qa are strongly affected since 
their relation is determined by redshift between 0.5 and 
1.8 (cf. Figure 1), where the systematic errors become 
larger, thus O m and Qa become limited by this type of 
systematic errors. An indication of systematic errors in 
Da from evolution might be found using an Q m - h plot 
(Figure 2b). 

In Figures 3a, b, and c we show similar constraints on 
fi m , w, and h from the SZ/X-ray method using the same 
fiducial CDM model, and 500 clusters with a random error 
of 4% in Da, as before (1, 2, and 3cr concentric ellipses, 
solid lines). The constraints form a banana-shaped region 
elongated in the Q m - w plane (Figure 3a). Systematic 
errors in the Da will not affect the results in the Q m — w 
plane, but simply shift the ellipses from random errors up 
or down along the h axis (Figures 3b and c, 3a solid el- 
lipses above and under the 1,2, and 3a concentric ellipses 
of random errors, solid lines). Again, with the assumed 
4% random and ±3% systematic errors, the h determina- 
tion is going to be limited by systematic errors. We also 
carried out simulations adding a systematic error with a 
gradient in redshift growing from 0% at z — to ±3% at 
z = 1 to the assumed 4% random error (Figures 3b and 
c, short dashed and dash dotted 3a lines). Again, we find 
that fi m , and w are strongly affected, their determination 
is limited by the assumed systematic errors. As before, 
the Hubble constant determination is not affected by this 
type of systematic errors. Also, an indication of system- 
atic errors in Da from evolution might be found using a 
w - h plot (Figure 3c). 

From Figure 2c, which shows constraints in the £Ia - h 
plane, we can conclude that with the assumed accuracy, 
the cosmological constant is not well constrained. From 
Figure 3a we can see that w is also poorly constrained by 
the SZ/X method. Clearly, whichever set of parameters 
is to be estimated, other cosmological measurements are 
needed to constrain these parameters further. 

One of the most promising other experiments to mea- 
sure cosmological parameters is based on CMB fluctua- 
tions. We estimated how well we can determine cosmo- 
logical parameters if we add constraints from CMB ex- 
periment, which we would expect to be particularly use- 
ful since the strongest dependency in this experiment is 
on the space curvature (the total average density in the 



Universe). As an illustration of the power of combining 
these techniques, we used the position of the first Doppler 
peak in the angular power spectrum of CMB fluctuations 
(Hanany et al. 2000; Mauskopf et al. 2000; Miller et al. 
1999), and the combination fl m h 2 — constant as observ- 
ables from primordial fluctuation studies (see for example 
Zaldarriaga, Spergel & Seljak 1997). 

The position of the first Doppler peak can be expressed 

as 

where k pea k and r(z*) are the first peak in k space and the 
effective distance to the last scattering surface at 



k( v (o) - V (z t )) 



(13) 



where K is the curvature, S is the sin, sinh, and identity 
function for models with negative, positive and flat space- 
time, and r] is the conformal time. In fc-space in a CDM 
model, 



kpeak ~ C\ + C 2 W m + C 3 W m - 

where w m = Q m h 2 and Wb 
are: c\ = 0.0112, c-i 



C4Wb + C 5 wl + C6W m Wb, (14) 

fif, ft 2 , and the coefficients 
0.0441, c 3 = -0.043, c 4 = -0.0496, 
c 5 = 2.65, and c 6 = 0.162 (White 1998). We used Hu & 
Sugiyama (1996)'s approximation for the redshift of the 
last scattering surface 

(15) 



z* = 1048[1 + 0.0012w-°- 7 ^][l + gf ] 



where 



and 



9i = 



0.0783w: 



1 + 39. 5w"; 763 ' 



.92 = 



0.560 



1 + 21.W- 8 



l ' 



(16) 



(17) 



In Figures 2a, b and c we over-plot the 3a error region 
based on the position of the first Doppler peak in the CMB 
power spectrum (long dashed lines). This region is roughly 
aligned along the line Q m + Q\ — constant because of the 
strong dependence of the position of the first Doppler peak 
on the total space curvature. For demonstration purposes, 



we choose £ v 



245 ± 10, where we assume a 3a error 



range of A£ = 10. As expected, constraints on cosmolog- 
ical parameters from SZ/X distance-redshift relation and 
the position of the first peak in the CMB fluctuations are 
highly complementary For our assumed set of clusters, 
uniformly sampled in z < 1, these two constraints are 
nearly orthogonal to each other in the Cl m - 51a plane (Fig- 
ure 2a). The constraints are also complementary in the 
h — fl m and h — Qa planes (Figures 2b and c) . Geometri- 
cally, we have narrow banana-shaped constraints from the 
SZ/X-ray method elongated in the Q m - tt\ plane and nar- 
row sheet-like constraints h = constant from the position 
of the first peak of the CMB fluctuation spectrum. The 
intercept of these constraints gives us stringent constraints 
on the following cosmological parameters: fi m , Q\, and h 
can be determined to ±0.08, ±0.1, and ±0.015 (3a, assum- 
ing 4% random error in Da)- Note, that while the effect 
on f2 m and fl\ from reasonable redshift independent sys- 
tematic errors are negligible, h would be dominated by 
systematic errors: a ±3% systematic error would result a 
±0.02 change in h. A systematic error which grows from 
0% to ±3% at z = 1 would cause an additional ±0.05, 
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±0.05, and ±0.005 error in Q m and and h. Con- 
straints from the SZ/X-ray method are also orthogonal to 
those from cluster evolution (compare constraints on the 
TO - S1a plane, our Figure 2a and Figure 1 of Holder et al. 
2001). In Figures 2a, b and c we also over-plot constraints 
from Q m h 2 = constant assuming a 10% error in its de- 
termination from CMB experiments (dash dot dot dotted 
lines), as suggested by studies based on the characteristics 
of the MAP experiment (Zaldarriaga et al. 1997). These 
constraints seem to be less useful in this parameter space 
than those from the position of the first Doppler peak. 

Constraints from the location of the first peak in the 
CMB fluctuations are not so useful when considering the 
alternative set of parameters w, Q, mi h. The I ~ 210 — 240 
constraint leads to surfaces which lie almost parallel to 
the ellipsoids derived from the SZ/X-ray method (com- 
pare our Figure 3a and Figure 4 of White 1998). In Fig- 
ures 3a, b, and c we over-plot constraints from fl m h 2 = 
constant again assuming a 10% error in its determination 
(dash dot dot dotted lines). Combining constraints from 
the SZ/X-ray method and those from Q m h 2 = constant 
(based on CMB fluctuation analysis), we obtain stringent 
constraints on w and h (< 0.2 and ±0.015, 3a). Note, 
that the effect on w from redshift independent systematic 
errors is negligible, and that, as before, h is dominated by 
systematic errors: a ±3% systematic error would result a 
±0.02 change in h. A systematic error which grows from 
0% to ±3% at z = 1 would cause an additional ±0.1 in 
w. From Figure 3b (Q TO — h plane), we see that a gradi- 
ent in the systematic error in redshift might be recognized 
using the £l m h 2 = constant constraint. Also, constraints 
from the shape of the power spectrum of CMB fluctua- 
tions (as will be achieved by MAP and Planck) are nearly 
orthogonal to those from the SZ/X-ray method (compare 
our Figure 3a and Figure 13 of Huterer and Turner 2000). 
Comparing our Figure 3a to Figures 8 and 9 of Haiman 
et al. (2000), we can see that constraints on the Q, m - 
w plane from distance function and cluster evolution are 
complementary. Also, constraints on w from the SZ/X- 
ray method combined with cluster abundance would put 
stringent constraints on w. 

Since the banana-shaped constraints from the SZ/X 
method in the parameter spaces defined by (f2 m , £!a, h) 
or (fl m , w, h) are nearly orthogonal to constraints from 
cluster evolution, a strategy for determining cosmological 
constants based only on clusters is likely to be successful. 
It seems to be possible to choose the parameters of the set 
of clusters used in this work such that constraints from the 
SZ/X ray method and cluster evolution can be made or- 
thogonal, and thus optimized for separating cosmological 
parameters. Clusters can provide a powerful, independent 
test for cosmological parameters. 

Note that our constraints in the fl m - w plane are curved 
towards the w axis (Figure 3a, solid and dotted lines), 
while Huterer & Turner (2000) 's constraints from SNe la 
(which is basically a constraints from distance - redshift 
function, cf. their Figure 23) show no sign of curvature. 
This is due to the fact that they used the Fisher matrix 
formalism, and we evaluated the % 2 statistic directly. Our 
results show that the likelihood function is strongly non- 
Gaussian on the fl m - w plane. 

As a second realization, we carried out simulations with 



seventy clusters (35 high redshift clusters, z > 0.5), a num- 
ber likely to be observed in the next few years, and esti- 
mated how well we can constrain cosmological parameters. 
We assumed the same fiducial cosmological model as be- 
fore: fl m — 0.3, f^A = 0.7, and h = 0.65. For this simu- 
lation, however, we assumed a random error of 7% in the 
angular diameter distance, which might be achievable in 
the near future (we assume lensing measurements will be 
used to eliminate errors from cluster shape). 

In Figures 4a, b, and c, we show the resulting 1, 2 and 
3cr constraints on the parameter space defined by (fl m , ^A, 
h), corresponding to A% 2 of 3.53, 8.02 and 14.2 appropri- 
ate for three fitted parameters (solid lines). Long dashed 
lines show the constraints from the position of the first 
Doppler peak. As before, we assumed l pea k = 245 ±10 (3ct 
range). We also over-plot constraints from Q m h 2 = con- 
stant assuming a 10% error in its determination from CMB 
experiments (dash dot dot dotted lines). From these Fig- 
ures we conclude that by using as few as 35 high-redshift 
(and 35 low-rcdshift) clusters, with a random error of 7% 
in Da we can constrain Q m , and h within ±0.2, ±0.2, 
and ±0.04 (3a errors). Note, that the effect of redshift in- 
dependent systematic errors in Da on f2 m and 17a , occur- 
ring in practice, is negligible, but a ±5% systematic error 
would result an additional ±0.035 error in h. Assuming 
that we know that our redshift independent systematic er- 
ror is less than 15% at the 3a level, with a 7% random error 
combined in quadrature, we would be able to determine h 
with an error of ±0.11. A systematic error which grows 
from 0% to ±5% at z = 1 would cause an additional ±0.1, 
±0.1 and ±0.01 error in J7 m , J7a, and h (Figure 4, short 
dashed and dash dotted lines). This means that as few 
as seventy clusters, with the errors likely to be achieved 
in the next few years, would be sufficient to exclude mod- 
els with zero cosmological constant with high significance. 
These results would be independent of the SNe la data, 
thus they would provide a robust check to the supernova 
results. 

In Figures 5a, b, and c, we show the 1, 2 and 3a con- 
straints on the parameter space defined by (fl„ l7 w, h), 
with the Q\ = 1 — il m constraint (solid lines). We over- 
plot constraints from £l m h 2 — constant assuming a 10% 
error as before (dash dot dot dotted lines). From Figure 5a 
we can see that we will be able to determine the equation 
of state parameter, w, with accuracy of 0.45 (3<r), even 
with reasonable systematic errors in Da, with additional 
information from CMB experiments, or SZ effect number 
counts. A systematic error which grows to ±5% at z = 
1 would cause an additional error of 0.2 in w (Figure 5, 
short dashed and dash dotted lines). 

5. CONCLUSION 

We have estimated the accuracy achievable in the de- 
termination of cosmological parameters using the SZ/X- 
ray method of distance determination. This method uses 
a sample of clusters to map the distance-redshift relation, 
which is a sensitive probe of cosmological parameters. The 
advantages of this well-known method are: (1) unlike other 
distance determination methods, it depends only on the 
geometry of the Universe, and its average densities; (2) it 
is a physical method, based on relatively simple gravita- 
tional virialization of clusters (as opposed to complicated 
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physics and chemistry involved in galaxy formation and su- 
pernova explosions); and (3) a large number of clusters is 
available for observation, and thus systematic effects can 
be reduced by using many clusters, or selecting clusters 
appropriately to reduce systematics. The necessary data 
should be available within the next few years. 

The SZ/X method, as other cosmological tests, can con- 
strain well only some combination of cosmological param- 
eters. We have shown that constraints on (f2 TO , ^a, h) 
from the Da{z) function measured in z = — 1 are nearly 
orthogonal to constraints from the first peak of the CMB 
fluctuations, and also to constraints from cluster evolu- 
tion. Constraints on (f2 TO , w, h) from the Da{z) function 
are complementary to those from cluster evolution (com- 
pare our Figure 2 to Figures 8 and 9 from Haiman et al. 
2000). In general, Da provides constraints on cosmological 
parameters similar to those from the SNe la method. Con- 
straints from cluster evolution (the source counts N$z{z) 
and Nx(z)) are similar to constraints from the full CMB 
fluctuation spectrum as will be measured by MAP and 
Planck. 

This result suggests that cosmological tests using only 
clusters of galaxies can be an important independent check 
on the values of cosmological parameters measured by 
other techniques. Cluster-based methods have systematic 
errors which are unrelated to those of other methods of 
measuring cosmological parameters (e.g., from CMB fluc- 
tuations, and SNe la). Furthermore, when the set of cos- 
mological parameters is further extended, with additional 
components of density or structure parameters, joint anal- 
yses using multiple cosmological tests will be essential to 
remove the parameter degeneracy exhibited by any one 
test. 

We demonstrated the effect of systematic errors on the 
determination of cosmological parameters. If random er- 
rors can be kept at a few percent level, systematic errors of 
similar magnitude will dominate the error in the Hubble 



constant. O m , and w are not affected by redshift- 
independent systematic errors. We also showed that a 
systematic error with a linear gradient in redshift will not 
affect Hubble constant determinations, but causes system- 
atic shifts in the estimates of f2 m , and w. 

We showed, that in the near future, even with only 35 
high redshift (and 35 low redshift) clusters, the SZ/X-ray 
method, combined with the position of the first peak in the 
power spectrum of CMB fluctuations, will provide enough 
accuracy to exclude 17a = models with high confidence 
even with constant and redshift dependent systematic er- 
rors (see Figure 4a). Also, the SZ/X-ray method, com- 
bined with the £l m h 2 — constant constraint, or cluster 
abundance, will allow us to determine the equation of state 
parameter, w, to within 0.45 (3<r; Figure 5a), even with 
usual (redshift-independent) systematic errors. 

This idealized discussion of the SZ/X-ray method and 
its errors in the determination of cosmological parameters 
could be improved by using the detailed characteristics of 
specific instruments and observing strategies. Departures 
from the assumptions made here could cause increases or 
decreases in the level of error in the derived parameters, 
with the largest changes likely for different redshift sam- 
plings. 

We conclude, that the determination of the angular di- 
ameter distance - redshift function using the SZ effect and 
X-ray thermal bremsstrahlung emission from clusters of 
galaxies can be used with confidence to constrain cosmo- 
logical parameters. In general, clusters of galaxies alone 
can be used to constrain cosmological parameters indepen- 
dently from other methods. Clusters lead to parameter 
limits which are competitive with using other techniques. 

Most of this work was done while SMM held a National 
Research Council Research Associateship at NASA God- 
dard Space Flight Center. We thank David Spergel and 
the anonymous referee for useful comments and sugges- 
tions. 
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Fig. 1. — Deviation of the angular diameter distance in Universes with Q m = 0.3 and Qa = 0.7, 0.6, and 0.3 relative to the angular distance 
function for a Universe with Q m = 0.3 and Qa = model, as a function of z (with h = 0.65). The solid, short dashed, and dash dotted lines 
represent percentage deviations with Qa = 0.7, 0.6, and 0.3 respectively. 
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Fig. 2a. — 1, 2 and 3cr constraints on cosmological parameters from simulated angular diameter distance measurements. We used a spatially 
flat CDM model as a fiducial model with Q m = 0.3, Qjy = 0.7, and h = 0.65. The figures show results of simulations using 500 clusters 
assuming a random error of 4% in the angular diameter distance (concentric solid ellipses). 2-dimensional (2D) projections of the 3D surface 
of 3cr constraints are shown using dotted lines. 3<r constraints assuming an additional ±3% systematic error are shown using solid ellipses 
above and below the random ellipses (Figure b and c). Constraints from a systematic error which grows from 0% to ±3% at z = 1 is shown 
using short dashed and dash dotted lines. Panels a, b, and c show two dimensional slices of the three dimensional parameter space defined by 
(f2 m , Ha, h) passing through the best-fit model as defined by the minimum of \ 2 ■ Constraints that might be derived from the location of the 
first Doppler peak, assuming £ pea k = 245 ± 10, are shown as long dashed lines, and define surfaces almost perpendicular to the error regions 
from Da(z) in the S7 m - h plane. Constraints from U m h 2 = constant, assuming a 10% error in its determination from CMB experiments, are 
shown using dash-dot-dot lines. 
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Fig. 2b. — 1, 2 and 3<x constraints on cosmological parameters from simulated angular diameter distance measurements. We used a spatially 
flat CDM model as a fiducial model with Q m = 0.3, SIa = 0.7, and h = 0.65. The figures show results of simulations using 500 clusters 
assuming a random error of 4% in the angular diameter distance (concentric solid ellipses). 2-dimensional (2D) projections of the 3D surface 
of 3<T constraints are shown using dotted lines. 3<r constraints assuming an additional ±3% systematic error are shown using solid ellipses 
above and below the random ellipses (Figure b and c). Constraints from a systematic error which grows from 0% to ±3% at z = 1 is shown 
using short dashed and dash dotted lines. Panels a, b, and c show two dimensional slices of the three dimensional parameter space defined by 
(n m , Ha, h) passing through the best-fit model as defined by the minimum of \ 2 ■ Constraints that might be derived from the location of the 
first Doppler peak, assuming £ pea k = 245 ± 10, are shown as long dashed lines, and define surfaces almost perpendicular to the error regions 
from Da(z) in the S7 m - h plane. Constraints from U m h 2 = constant, assuming a 10% error in its determination from CMB experiments, are 
shown using dash-dot-dot lines. 
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Fig. 2c. — 1, 2 and 3<r constraints on cosmological parameters from simulated angular diameter distance measurements. We used a spatially 
flat CDM model as a fiducial model with Q m = 0.3, Qjy = 0.7, and h = 0.65. The figures show results of simulations using 500 clusters 
assuming a random error of 4% in the angular diameter distance (concentric solid ellipses). 2-dimcnsional (2D) projections of the 3D surface 
of 3cr constraints are shown using dotted lines. 3<r constraints assuming an additional ±3% systematic error are shown using solid ellipses 
above and below the random ellipses (Figure b and c). Constraints from a systematic error which grows from 0% to ±3% at z = 1 is shown 
using short dashed and dash dotted lines. Panels a, b, and c show two dimensional slices of the three dimensional parameter space defined by 
(n m , Ha, h) passing through the best-fit model as defined by the minimum of \ 2 ■ Constraints that might be derived from the location of the 
first Doppler peak, assuming £ pea k = 245 ± 10, are shown as long dashed lines, and define surfaces almost perpendicular to the error regions 
from Da(z) in the S7 m - h plane. Constraints from U m h 2 = constant, assuming a 10% error in its determination from CMB experiments, are 
shown using dash-dot-dot lines. 
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Fig. 3a. — 1, 2 and 3<x constraints on cosmological parameters from simulated angular diameter distance measurements. Contours and 
lines have the same meaning as in Figure 2. The figures show results of simulations using 500 clusters assuming an error of 4% in the angular 
diameter distance (concentric ellipses, solid lines). 3<r constraints assuming an additional ±3% systematic error are shown using solid ellipses 
above and below the random ellipses (panels b and c). Panels a, b, and c show two dimensional slices of the three dimensional parameter 
space defined by (f2 m , w, and h) and passing through the best-fit model as in Figure 2. Constraints from fl m h 2 = constant, assuming a 10% 
error in its determination from CMB experiments, are shown using dash-dot-dot lines. 
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Fig. 3b. — 1, 2 and 3cr constraints on cosmological parameters from simulated angular diameter distance measurements. Contours and 
lines have the same meaning as in Figure 2. The figures show results of simulations using 500 clusters assuming an error of 4% in the angular 
diameter distance (concentric ellipses, solid lines). 3<r constraints assuming an additional ±3% systematic error are shown using solid ellipses 
above and below the random ellipses (panels b and c). Panels a, b, and c show two dimensional slices of the three dimensional parameter 
space defined by (f2 m , w, and h) and passing through the best-fit model as in Figure 2. Constraints from fl m h 2 = constant, assuming a 10% 
error in its determination from CMB experiments, are shown using dash-dot-dot lines. 
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Fig. 3c. — 1, 2 and 3cr constraints on cosmological parameters from simulated angular diameter distance measurements. Contours and 
lines have the same meaning as in Figure 2. The figures show results of simulations using 500 clusters assuming an error of 4% in the angular 
diameter distance (concentric ellipses, solid lines). 3<r constraints assuming an additional ±3% systematic error are shown using solid ellipses 
above and below the random ellipses (panels b and c). Panels a, b, and c show two dimensional slices of the three dimensional parameter 
space defined by (f2 m , w, and h) and passing through the best-fit model as in Figure 2. Constraints from fl m h 2 = constant, assuming a 10% 
error in its determination from CMB experiments, are shown using dash-dot-dot lines. 
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Fig. 4a. — Same as Figure 2, but with only 70 clusters (with 35 at z > 0.5) assuming a 7% random error in the angular diameter distances. 
We omit the redshift-independent systematic error ellipses, since they simply correspond to the same shifts on the h axis seen in Fig. 2, but 
retain the redshift-dependent systematic error ellipses corresponding to a 5% gradient in the error in Dj\ . 
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Fig. 4b. — Same as Figure 2, but with only 70 clusters (with 35 at z > 0.5) assuming a 7% random error in the angular diameter distances. 
We omit the redshift-independent systematic error ellipses, since they simply correspond to the same shifts on the h axis seen in Fig. 2, but 
retain the redshift-dependent systematic error ellipses corresponding to a 5% gradient in the error in Dj\ . 
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Fig. 4c. — Same as Figure 2, but with only 70 clusters (with 35 at z > 0.5) assuming a 7% random error in the angular diameter distances. 
We omit the redshift-independent systematic error ellipses, since they simply correspond to the same shifts on the h axis seen in Fig. 2, but 
retain the redshift-dependent systematic error ellipses corresponding to a 5% gradient in the error in D^. 
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Fig. 5a. — Same as Figure 3, but with only 70 clusters (with 35 at z > 0.5) assuming a 7% random error in the angular diameter distances. 
We omit the redshift-independent systematic error ellipses, since they simply correspond to the same shifts on the h axis seen in Fig. 2, but 
retain the redshift-dependent systematic error ellipses corresponding to a 5% gradient in the error in Dj\ . 
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Fig. 5b. — Same as Figure 3, but with only 70 clusters (with 35 at z > 0.5) assuming a 7% random error in the angular diameter distances. 
We omit the redshift-independent systematic error ellipses, since they simply correspond to the same shifts on the h axis seen in Fig. 2, but 
retain the redshift-dependent systematic error ellipses corresponding to a 5% gradient in the error in Dj\ . 
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Fig. 5c. — Same as Figure 3, but with only 70 clusters (with 35 at z > 0.5) assuming a 7% random error in the angular diameter distances. 
We omit the redshift-independent systematic error ellipses, since they simply correspond to the same shifts on the h axis seen in Fig. 2, but 
retain the redshift-dependent systematic error ellipses corresponding to a 5% gradient in the error in D^. 



